rm(list = ls(all = TRUE)) #clear workspace

library(foreign)
library(survey)
options(survey.lonely.psu = "adjust")

setwd("C:/Users/martin.romero/Dropbox/ENVIPE/2011/Base_de_datos")
#setwd("~/Dropbox/ENVIPE/2011/Base_de_datos")

################################################
#Pregunta al entrevistado
################################################

data <- read.dbf("tper_vic.dbf")

#Ponderador y diseno
data$FAC_ELE <- as.numeric(as.character(data$FAC_ELE))
design.per <- svydesign(id=~UPM,strata=~EST, weights=~FAC_ELE, data=data)

#Lo secuestraron para exigir dinero o bienes:
# 1 Si
# 2 No
# 8 No responde
# 9 No sabe
table(data$AP8_3_11) #77 casos
foo <- svymean(~AP8_3_11, design=design.per, na.rm=TRUE) #0.12%
print(foo)
confint(foo)
foo <- svytotal(~AP8_3_11, design=design.per, na.rm=TRUE) #88,784
print(foo)
confint(foo)

#Me puede decir cuantas veces?
table(as.numeric(as.character(data$AP8_4_11)))
sum(as.numeric(as.character(data$AP8_4_11)), na.rm=TRUE) #86
foo <-svytotal(~as.numeric(as.character(AP8_4_11)), design=design.per, na.rm=TRUE) #114,958
print(foo)
confint(foo)

################################################
#Modulo de victimizacion
################################################

modulo <- read.dbf("tmod_vic.dbf")

#7.1 �Me podr�a decir si sus secuestradores.  					
#1  exigieron rescate a familiares por su liberaci�n?
#2	"lo obligaron a retirar dinero de un cajero, entregar joyas,celular u otras cosas?"
#8	No responde
#9	No sabe
#7	No especificado 
#b	blanco
table(modulo$BP7_1)
round(prop.table(table(modulo$BP7_1)),3)*100

#7.2 �Cu�nto tiempo lo tuvieron secuestrado?  					
#1  Menos de 24 horas
#2	De 1 a 3 d�as
#3	De 4 a 29 d�as
#4	De 1 a 3 meses
#5	M�s de 3 meses
#8	No responde
#9	No sabe
#7	No especificado 
#b	blanco
table(modulo$BP7_2)
round(prop.table(table(modulo$BP7_2)),3)*100













